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Abstract 

We study the state-counting problem that arises in the SU{2) black hole 
entropy calculation in loop quantum gravity. More precisely, we compute the 
leading term and the logarithmic correction of both the spherically symmetric 
and the distorted SU{2) black holes. Contrary to what has been done in 
previous works, we have to take into account "quantum corrections" in our 
framework in the sense that the level k of the Chern-Simons theory which 
describes the black hole is finite and not sent to infinity. Therefore, the new 
results presented here allow for the computation of the entropy in models 
where the quantum group corrections are important. 



Introduction 

The black hole entropy calculation in the framework of loop quantum gravity [1] is 
based on the effective description of the quantum gravitational degrees of freedom 
at the black hole horizon obtained from a suitable quantization of the classical phase 
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space describing isolated horizons (see p] and references therein). In these models 
the degrees of freedom at the horizon are described by Chern-Simons theories with 
SU{2) (or f/(l)) structure groups [SI SI El El E]- The simplest models are those 
where spherical symmetry is imposed already at the classical level. In this case it is 
natural (although not necessary) to consider SU{2) (or U{1)) Chern-Simons theory 
with a level that scales with the macroscopic classical area /c oc a^. This makes 
the state-counting (necessary for the computation of the entropy) a combinatorial 
problem which can be entirely formulated in terms of the representation theory of 
the classical group SU(2) (or U{1)): for practical purposes one can take k = oo from 
the starting point [HI El [IDl [II] • 

However, the perspective considered above can be completely changed if one 
studies the models in the recently introduced SU{2) Chern-Simons formulation 
[SI [SI [Zj- The necessity of an SU{2) gauge invariant formulation comes from the 
requirement that the isolated horizon quantum constraints be consistently imposed 
in the quantum theory (in [5j it is shown how the U{1) treatment leads to an artifi- 
cially larger entropy due to the fact that some of the second class constraints arising 
from the SU{2)-to-U{l) gauge fixing can only be imposed weakly). However, the 
SU{2) formulation is not unique as there is a one parameter family of classically 
equivalent SU{2) connections parametrizations of the horizon degrees of freedom. 
More precisely, in the passage from Palatini-like variables to connection variables 
that is necessary for the description of the horizon degrees of freedom in terms of 
Chern-Simons theory (central for the quantization), an ambiguity parameter arises 
[HI E]- This is completely analogous to the situation in the bulk where the Immirzi 
parameter reflects an ambiguity in the choice of SU (2) variables in the passage from 
Palatini variables to Ashtekar-Barbero connections (central for the quantization in 
the loop quantum gravity approach). In the case of the parametrization of the iso- 
lated horizon degrees of freedom, this ambiguity can be encoded in the value of 
the Chern-Simons level fc, which, in addition to the Immirzi parameter, becomes an 
independent free parameter of the classical formulation of the isolated horizon-bulk 
system. 

Therefore, it is no longer natural (nor necessary) to take k oc Qh- On the con- 
trary, it seems more natural to exploit the existence of this ambiguity by letting k 
be arbitrary 0. In this way the SU{2) classical representation theory involved in pre- 
vious calculations should be replaced by the representation theory of the quantum 
group Uq{su{2)) with q a non-trivial root of unity. Thus quantum group corrections 
become central for the state-counting problem. In this paper we study the finite k 
counting problem by means of simple asymptotic methods. The powerful methods 
that have been developed for the resolution of the counting problem in the k = oo 
[TU| [TT] are perhaps generalizable to the finite k case. Here we follow a less rigorous 
and more physical approach. The formulation is partly inspired from a combination 

^ It is a good thing that the effective treatment contains a free parameter arising at the boundary 
from exactly the analogous reason as the Immirzi parameter in the bulk parametrization of the 
phase space. This keeps open the possibility that dynamical considerations could lead to cancelation 
of both ambiguities producing Immirzi parameter independent predictions. 
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of ideas stemming from different calculations in the literature [9l [121 US] • 

In the first section, we review some basic facts concerning the quantization of 
SU{2) Chern-Simons theory whose physical states are built from the representation 
theory of the quantum group Uq{su{2)) when g is a root of unity. For that reason, 
we recall some properties of the representation theory of Uq{su{2)), which allows 
us to compute the dimension of the Chern-Simons theory Hilbert space Tics when 
the space is a punctured two-sphere. In the second Section, we give a new integral 
formulation of the dimension of Tics which appear to be much more convenient to 
compute black hole entropy. In the last Section, we compute the leading term of 
the S'[/(2)-black hole entropy and its logarithmic corrections first for the spherically 
symmetric black hole and then for the distorted black hole. We adapt the techniques 
used in [TU] and firstly introduced in [Sj to compute the entropy of a black hole. We 
are not going into the mathematical details of these techniques which has been 
very well exposed in [lOJ and which are in fact very well-known in the domain of 
probabilities and used to understand some properties of random walks. We recover 
that in the spherically symmetric and distorted black holes, the leading term of the 
entropy is proportional to the area and the first corrections are still logarithmic: 
S{a) ~ aa -f- (3\oga. In the spherically symmetric case, a depends on the level 
k and reaches the value obtained in previous calculations when k goes to infinity; 
concerning (3 it is independent of k and is given by the value —3/2 as expected. In 
the distorted case, a grows logarithmically with k and f3 is fixed to the value —3. 
We finish with a discussion. 

1. The Chern-Simons Hilbert Space 

The Chern-Simons theory associated to the group SU{2) is a gauge theory on a 
three dimensional manifold M governed by the action 



where k is called the level of the action, A is the local S'f/(2)-connection field and 
(•, •) is a notation for the 5u(2) Killing form. The Chern-Simons theory became 
really important when first it was shown to be closely related to three dimensional 
gravity [HI |15] and above all when Witten showed [16] its amazing relation to 
manifold and knots invariants. Indeed, the Chern-Simons path integral is a manifold 
invariant whereas the mean values of quantum observables naturally lead to Jones 
polynomials. For all these reasons, Chern-Simons theory has been the center of a 
lot of interests and its quantization is now very well-known when the gauge group 
is compact, and in particular when the gauge group is SU{2). 

The covariant (path integral) and canonical quantizations offer the two main 
strategies to quantize the Chern-Simons theory. These approaches are complemen- 
tary: the covariant quantization leads easily to the fact that the level k must be an 
integer when the gauge group is compact [16]; the canonical quantization leads to 
a precise description of the Hilbert space when the gauge group is compact but not 




(1) 
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only (see [T7] for an introduction of the combinatorial quantization for example). 
Both quantizations are necessary to understand how the mean value of Wilson loops 
observables are related to knots polynomials like the (colored) Jones polynomial or 
its generalizations. 

Here we are exclusively interested in the description of the Hilbert space of 
Chern-Simons theory when the space is a two-sphere punctured with a number p 
of particles. At the classical level, each puncture, labelled by £ G comes with 

an unitary irreducible representation ji of the gauge group SU{2). At the quantum 
level, one shows that the classical group gauge symmetry is replaced by a quantum 
group symmetry and the Hilbert space is constructed from the representation theory 
of the quantum group Uq{su{2)) where q = exp{i7i/{k + 2)) is necessarily a root 
of unity. An immediate consequence is that the SU{2) representations labeling the 
classical punctures become Ug{su{2)) representations which concretely implies a cut- 
off on the punctures' representations which cannot be higher than k/2. Then, the 
associated Hilbert space is the vector space 



of invariant tensors in the tensor product <^eje of Uq{su{2)) representations endowed 
with a Hilbert structure defined, as in the classical case, from the Haar measure on 
the quantum group. However, we will not be interested in the Hilbert structure of 
Hk{ji, ■ ■ ■ in the rest of the paper but rather in its vector space structure and 
more precisely in its dimension. Indeed, the computation of the SU{2) black hole 
entropy in Loop Quantum Gravity needs to be done precisely the calculation of the 
dimension of the previous vector space. 

The calculation of the Hilbert space dimension makes use of the Verlinde coeffi- 
cients. In order to introduce these coefficients, we start by recalling some basic facts 
concerning the representation theory of Uq{su{2)). 

1.1. Basics of the representation theory of Uq{su(2)) 

This Section is devoted to recall some basic results on the quantum group Uq{su{2)) 
we will need in the sequel. We are not going to give a precise definition of this 
quantum group and a complete description of its properties. Furthermore we will 
be interested only on some aspects concerning its representations theory and its 
recoupling theory. 

The (standard) unitary irreducible representations of Uq{su{2)) are labelled by 
integers j < k/2. The dimension dj of the j-representation is the same as in the 
classical theory and then we have dj = 2j + 1. Given an element ^ G Uq{su{2)), 
its representation vrj(^) is an endomorphism of the vector space Vj. Many formulae 
coming from the representation theory of Uq{su{2)) are expressed in terms of q- 
numbers [x] defined for any complex number x by the relation: 



HkUi, ■■■ ,jp) = Irivi^eje) 



(2) 



[x] 



— q 
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Invariant f/g(s-u(2))-tensors are defined, by analogy with the classical situation, as 
tensors which are invariant under the adjoint action. Note however that the adjoint 
action is deformed compared to the classical case and makes use of the antipode 
instead of the inverse. Among the invariant tensors, the 3-valent ones 



are particularly interesting because all invariant tensors decompose into 3-valent 
intertwiners. Three-valent intertwiners are represented as usual by a vertex be- 
tween three lines colored by the representations je as illustrated in the figure (P). 
Contrary to what happens in the classical case where one can make a certain choice 



Figure 1: Pictorial representation of the 3-valent intertwiner i(ji,j2;j3) and its 
adjoint operator i^s; ji, j2)- 

of normalization such that the matrix elements of J2; Js) are invariant under the 
permutations of (ji, ^2)5 the order between the representations in t(ji,j2,i3) does 
matter because of the presence of a non-trivial braiding in the quantum case. 

In the sequel, we fix the normalization of t{ji,j2', js) such that it satisfies the 
following fusion rule: 



where Ij^i^ja is the identity map in the tensor product of the two representations ji 
and j2 and ^(js; ji, ^2) : Vj3 Vji ® Vj^ is the adjoint of i-iJi, j2i js)- This relation 
implies that the so-called ^-graph is normalized to one 



Trii®j2(i(ji, j2; j3)-i(j3; ji, j2)) = Tr^sl^Oa; ji, j2)-i(ji, j2; j3)) = >"(ji,j2,i3) , (4) 



We made used of the notation I'O'i, J2, J3) £ {0, 1} which is one only when (ji, J2, J3) 
satisfy the triangular inequalities, otherwise it vanishes. These identities are graph- 
ically represented in the figure ([2]). 

Of particular interest for the quantization of Chern-Simons theory is the fact that 
Uq{su{2)) is quasi-triangular and therefore admits an universal R-matrix which is 
at the origin of the braiding properties associated to the quantum groups. Without 
entering too much into the details, let us recall that R G Uq{su{2)) (g) Uq{su{2)) 






(3) 



which is equivalent to 
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Figure 2: Illustration of the normalization properties: the three relations are in fact 
equivalent to the condition that the ^^-graph is normalized to one. 

satisfies in particular the so-called quantum Yang-Baxter equation and other defining 
properties that one can find in [18] for example. 

The evaluation of the R-matrix in the tensor product of representations ji ® j2 
is denoted -Rjija = (^ii ® '^j2){R) ^'^^ defines a braiding operator from Vj^ ® Vj^ to 
the opposite tensor product Vj^ ® Vj^^. It is useful to represent the R-matrix as in 
the picture ([S]): if i? is represented by an under-crossing (the up-line undercrosses 
the down-line) then its inverse is represented by an over-crossing (the up-line 
overcrosses the down-line). It is clear from this representation that the product of 
R by its inverse is the identity because the braiding has been unknoted. 

Jl 32 Jl y ^2 



32 Jl 32 Jl 

R 

Figure 3: Pictorial representation of the R-matrix and its inverse R^^. Both R- 
matrices are evaluated in ji ® j2- 



1.2 From Verlinde coefRcients to the Hilbert space dimension 

Now, we have all the ingredients to construct the Verlinde coefficients. These coef- 
ficients appeared first [19] in the context of conformal field theory and then it has 
been realized that they have a very simple algebraic interpretation in the context of 
quantum groups. Here we will give only their algebraic definition and some of their 
properties which are important in the calculation of the dimension of the Hilbert 
space Hk{ji, ■■■ ,je). 

Given two unitary irreducible representations ji and j2, one defines the Verlinde 
coefficient Sj-^^j^ = Sj^j^ as the real number determined by the trace on Vj^ ® Vj^ of 
the operator R^ up to a normalization factor: 

^jij2 — ^ Tr (i^jjj^-Rjijj) (6) 
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where 

Z = J-^sm(-^) (7) 

is in fact the partition function of the SU (2) Chern-Simons theory on the 3-sphere 
S^. It will be useful in the sequel to use the "un-normalized" Verlinde coefficient 
Sjij2 = Tr (i?j2ii-Riii2) ^^e choice of the normalization factor will appear clear 
soon. Note that Sj-^^j^ is the evaluation on the Hopf-link embedded into the 3-sphere. 
The Hopf-link is represented in the figure (jl]). 




Figure 4: Representation of the Hopf-link. The evaluation of the associated quantum 
spin-network colored with the representations ji and j2 gives the un-normalized 
Verlinde coefficient Sj^^j^. 

The explicit expression of the R-matrix implies that 

sinf^^^^ 



^jlj2 ~ [d'jld'j2] - 



These coefficients satisfy many interesting properties which are important to 
compute the dimension of the physical Hilbert space iffcOi, ■ ■ ■ ,je) presented above. 
The properties we will need are given below: 

The normalization relation: ^jij3^j3j2 = ^2^' 

33 

The fusion relation: Sj^^j^Sj^j.^ = [dj^]''^^Y{ji, j2, i) Sj.^g] (10) 

e 

n—l n— 1 

The recursive relation: JJ Sj^j^ = ^ 5ei,o JJ yijiJiJi+i)SE„j„ill) 



The definition of the normalized Verlinde coefficient becomes clear from the normal- 
ization relation. The recursive relation is a generalization of the fusion relation to 
any number of unitary irreducible representations j = (ji, ■ ■ ■ ,jp)- The fusion and 
recursive relation are really easy to prove using the graphical representations of the 
Verlinde coefficients. The proof of the fusion relation is given in the picture ([5]); the 
proof of the recursive relation is done along exactly the same lines. 

Verlinde coefficients and their properties are particularly interesting to obtain 
useful formulae for the dimension of the Hilbert space H^^ji, ■ ■ ■ ,jp). Indeed, the 
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Figure 5: Pictorial proof of the fusion relation. We start with the graph on the left. 
The two arrows are identities: the first one is obtained applying the decomposition 
of the identity "along the vertical dashed hne" ; the second one is obtained applying 
the decomposition of the identity "along the horizontal dashed line". Both lead to 
equivalent expressions and the equality between the evaluations of the graphs on the 
right is exactly the fusion relation. We made used of the identities represented in 
the picture ([2]). 

dimension A^fc(j) = dim(i^fc(ji, ■ ■ ■ ,jp)) is 

p 

ii,---,ep 1=1 

and can be expressed in terms of Verlinde coefficients using the recursive relation 
(|TT|) combined with the normalization relation Some trivial calculations lead to 
the expression: 

which reduces, after using the explicit formula of Verlinde coefficients ([H]), to the 
following well-known formula: 

Nkij) = T^y(sin(-^))^-^fTsin(^). (14) 

I 1=1 

2. Equivalent formulae for the Hilbert space dimension 

The expression ([T3j) for the dimension of the SU{2) Chern-Simons Hilbert space is 
not very useful to compute the entropy of a Black Hole. We propose here to give 
equivalent more interesting formulae. 

2.1. Chern-Simons Hilbert space and random walk 

The fact that the coefficient Nk{j) are closely related to random walks have been 
noted and investigated in [13] in the classical case, namely when k becomes infinite. 
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Here we show that even in the quantum case (i.e. for a finite k) this hnk between 
Chern-Simons and random walk still exists and appears very interesting for the 
calculation of the entropy. Some of our formulae have been derived in [12] where 
however only the classical case (infinite k) has been studied at the end. For obvious 
reasons of notations, we will consider Nk{d) = Nk-2{3) in the sequel. We will also 
make use of the notations di = dj. = 2ji + 1 and d = {di, - ■ ■ , dp). 

This Section is devoted to propose a random walk interpretation of the dimension 
Nk{d). For that purpose, we proceed in four steps. 

a) A^fe(d) as a difference of i?fc(d)-type functions 

To do so, we first use the identities sin^ 9 = 1 — cos^ 9 and cos 9 = sin(26') / (2 sin^ 9) 
in the formula 

^,.) . 5:,.(^),n555^ (15, 

to write it as the difference 

iV,(d) = Bk{d) - i5,(d+) (16) 

where the function Bk which depends on a family of representations d = {di, ■ ■ ■ ,dp) 
or d+ = {di, ■ ■ ■ ,dp,2, 2) reads 

2 sin(^rf£) , , 

The product runs over i G [0,p] or £ G [0,p + 2] depending whether we are consid- 
ering -Bfc(d) or Bk{d+). Note that d+ is the union of the family of dimensions d 
with two more equal elements corresponding to the dimension of the fundamental 
representation di/2 = 2. 

b) Combinatorial expression of Bk{d) 

Now we concentrate on the function Bk{d). In particular, we want to exhibit the 
fact, as in the classical case, that Bk{d) admits a random walks interpretation. To 
show this is indeed the case, we replace each term of the product in (1171) by the 
following expression: 

^'""^ k = e^f (^.-1) V e-2^ir-^ (18) 
sin(^) ' ^ ^ 

As a result, the function ffTTl) can be rexpressed as follows: 

^^(d) = ^EriE^''^^''"'"'"'^ = ^E E e'"^^^'^-'""^ (19) 

d=0 1=1 ni=0 d=0 {m,--- ,np} 
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where we introduced the notations Ap = ~ 1) ^ — S?=i Note that 

the second sum runs over famihes of integers {ni, ■ ■ ■ , Up} such that each component 
rii G [0, di — 1]. Permuting the two sums and summing over the variable d lead to: 

{r!,i,---,np} {ni,---,np} 

To go further into the calculation, we distinguish the case where is odd from the 
case where is even. 

c) Af odd implies that A^fc(d) = 

The first case, A^ odd, is simpler. Indeed, in that case, 1 — e™^'^ = 2 and then the 
function Bk{d) reduces to the form: 

{ni,--- ,np} 

From the beginning, we know that -Bfc(d) is a real- valued function and therefore it 
equals its real part, i.e. Bk{d) = 7le{Bk{d)) where 7le{z) denotes the real part of 
z G C. Moreover, for any value of 6 (different from 0[27r]), the following equality 
holds: 

1 1 i e 

= - + -cotan- . (22) 

1 - e*^ 2 2 2 ^ ' 

As a consequence, the expression of the function Bk{d) simplifies drastically and 
reduces to: 

^^(d) = ^^ J2 1 = Iflde. (23) 

{ni ,np} £=1 

Therefore, i?fc(d_|_) = 4Bk{d) and then the dimension of the Hilbert space f lT^ 
vanishes in that case. The meaning of this result is simple: there is no invariant 
tensor in the tensor product ^ije when Ap = J^ei^^i ~ 1) is odd. As an example, let 
us consider the case where all the spins equal 1/2: Ap = p odd means that there is 
an odd number of spins; as expected there is no trivial representation in the tensor 
product of an odd number of 1/2 representations. 

d) Ai even: random- walk interpretation of A^A:(d) = 

The second case, Ap even, is far more interesting. In that situation, we would naively 
say that -Bfc(d) vanishes because all the terms in the numerator of the formula ( !20l) 
are 1 — e*'^^*' = 0. But a more careful analysis shows that the denominator can also 
lead to a singularity. As a result, the non- vanishing contributions to the sum ( l20l) 
are those where both the numerator and denominator vanish. For this to happen, 
there must exist s G Z such that Ap — 2A^ = 2sk. Therefore, we have: 

Bk{d) = 2 ^Ap-2N[2k] (24) 

{jii,--- ,np} 
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where 6n[2k] takes the value one if there exists an integer s such that n = 2ks, 
otherwise it is null. Here comes the random walk interpretation of the dimension 
of the Chern-Simons Hilbert space. We proceed to a changing of variables: instead 
of summing over non-negative integers rii G [0,di — 1], we sum over half-integers 
rrii = rii - ^ G (with mj+i = + 1). Then, the formula (Elj) becomes 

-Sfe(d) = 2 ^ 6mj_+-+mj,[k] ■ (25) 
{mi,--- ,mp} 

A similar formula has been found in [12] and its classical counterpart has been 
established and studied in [I3]. As a result, the Bk function appears to be the 
number of ways to start from the origin at the Z axe and come back at a point 
0[k] after p steps mj, each step being bounded as follows G [^-^, ^-^]. These 
functions have been deeply and precisely studied in the domain of random walks. 
And this analogy was used as a central tool in [13] to obtain asymptotics behavior 
of some entropy. In order to be a bit more explicit, we introduce the variable 
r = [Ap/ {2k)] where [x] is the floor function and then: 

r 

■■+mp—qk • 

(26) 

{mi,---,mp} q=-r 

Let us recall now that we are interested in the number of states and not in the 
function itself. Using previous formulae, we have for A^^ the expression: 



/ 1 

A^fe(d) = 2 ^ ^ (5mi+---+mp-gfe - - E ^mx 



{mi,---,mp} \q=-r a,b&{-\,\}<l 



^ ^ ^ ^ - \-m.p+a+h—qk 



where s = [(Ap + 2)/ {2k)]. It is clear that s belongs to the set {r, r + 1} and to avoid 
complications we assume that r = s. The case s = r + 1 would introduce extra terms 
which are not important at all for what we want to do. In that case, the previous 
formula simplifies and after summing over the variables a and b one obtains: 

~ 1 1 

A^fc(d) = {6mi+-+mp-qk — ■^Smi + ---+mp-qk+l — -^Smi+-+mp~qk-l) ■ (27) 

{mi,--- ,mp} q=-r 

This expression generalizes the one obtained in the classical case (which corresponds 
in fact to r = in our notations). It is useful to study the asymptotic behavior of 
the number of states and also to study the effect of a finite k. 



2.2. Integral formula 

Very often, one identifies the number of states A^^ to the dimension of the invariant 
tensors space in the tensor product ^eje between representations of the classical 
group SU{2). This is only true when the ratio r = which also coincides with 
the classical limit k goes to infinity. In that case, A^^ = A^oo is expressed as an 
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integral over SU{2) conjugacy classes, or equivalently over an angle 6. We proposed 
to generalize this integral formula to the case where r 7^ 0. 
For that purpose, we start with the relation: 







27r 

defined for any integer a G Z. It easily leads to the relation: 

> Sm,+-+mj,+a d^COs(a^)|| 

{mi,---,mp} <:=1 ^ 



(29) 



where the sum runs over G [^-3^, ^-^]. Using this last identity and after some 
trivial calculations, one shows that the number of states is given by the integral: 



N,{d) = - f^dOij^ cos(^#)) sin2(^) n . (30) 

One more simplification occurs due to the trigonometric identity 

cos{eqk) = 1 + 2 ^ cos(gA;^) = ^^(^ +^2)^^) _ ^3^^ 

sm TT 

q=—r q=l ^ 

As a result, the number of states Nk takes the form: 

'■^'^ sin((r + l)ke) sin(rf4) 



This formula generalizes, as announced in the introduction of that section, the clas- 
sical one. We see, as expected, that A^fc(d) coincides with the classical formula when 
r = 0, i.e. when Ap < 2k. This particular case can be recovered from different argu- 
ments: if Ap < 2k, then each representations in the tensor product <S>ije has a spin 
s < k/2 and therefore one never sees the effect of the cut-off k. When the condition 
r = is not satisfied the integral formula defining the number of states differs from 
the classical one by a different integration measure on the SU{2) conjugacy class. 
The presence of this new measure might have an effect on the black hole entropy. 



2.3. An example: all spins equal 1/2 

To get an intuition of previous formulae, we consider a particular example: we 
assume that dj^ — 2 for all punctures i G {1, • • • ,p}, i.e. all the spins equal 1/2. 
Furthermore, we assume that Ap = p is even. 
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a) Classical case: r = |13|, 

To start with, we also consider first tlie case Ap < 2k, i.e. r = in the previous 
notations. This case has been studied deeply in [T3] . 

From the random walk expression of the number of states (1271) . one obtains that 
Nk{2,--- ,2) = N2{p) is given by: 

= (^rn,+-+m, " S^,+...+m,+l) = (^^^^ " l) ' ^^^^ 

{ TTl 1 , ■ ■ ■ , TTLp J- 

We have omitted to mention k because Nk does not in fact depend on k when r = 0. 
The last equality is a result of a trivial combinatorial analysis: given an integer 
a < p/2, the number of ways to have rii + ■ ■ ■ + Up = a where G {—1/2, 1/2} is 
given by the number of ways to choose {p/2 + a) elements from a set of p elements 
which is precisely given by the binomial coefficient (p/2+a)' From the expression 
of the binomial coefficients in terms of factorials, one ends up with the following 
formula: 

'^^l^' - ^2jmV '''' 
We obtain an exact combinatoric formula for the number of states in that particular 
case. The asymptotic of N2{p) is therefore straightforward to obtain from the Stirling 
formula which states that: 

p\ ~ \/2'np[j)l ef for large values of p. (35) 

Using this very well-known result, one shows the following asymptotic behavior: 



Mv) ~ \ -P''^'2P. (36) 



This formula coincides with the one found in [12] from different arguments. In 
particular, we recover the same leading order and the same sub-leading corrections 
to the "entropy": 

S,/2ip) ^ \0gN2ip) = 2^ - hogp + 0(1). (37) 
b) Quantum corrections: r > 0. 

Let us now relax the condition that Ap < 2k, i.e. r = [p/{2k)] is now a non-zero 
integer. In that case, it is a bit more involved to obtain a combinatoric expression 
for the number of states but, using similar arguments as previously, one can show 
that: 



^2(p) 5Z 2 (j9/2 - gA; - l) 2 (p/2 - gA; + l) 



(3^ 
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Note that we still omit to mention explicitly the dependence in k even if now N2{p) 
does depend on k. To go further, we separate the g = contribution from the others 
in the sum and, using trivial symmetries properties of binomial coefficients, we get: 



P ] _ ( P 
p/2j {p/2 - 1 



9=1 



p/2 - qk) \p/2 -qk-l) \p/2 - gA; + 1 ' ' 



Thus, we obtain a correction to the previous classical case fl33p and each term (for 
a given value of g G [l,r]) in the remaining sum is a finite linear combination of 
binomial coefficients which reduces to the following form after some calculations: 

1-p- 2g2A;2 f p ^ ^^^^ 



{p/2 + 1)2 - g2p \^p/2 - qk 

As a result, the asymptotic behavior of each term in the previous sum is governed 
by the asymptotic behavior of the binomial coefficient (p/2-gA:) ■ interested in 

the asymptotic for a large value of p but also a large value of k such that these two 
numbers have the same scaling namely p/k remains constant. Indeed, in the black 
hole context, both p and k are proportional to the area of the horizon which tends to 
infinity (in unit of Planck area). Thus, the leading term in the asymptotic expansion 
of (HOl) is given by the leading term of the expansion of binomial coefficient of the 
form (^p) for large values of p where a E [0, l/2[. From Stirling formula (!35ll . one 
shows that 

^ g{a)-P with ^(a) = a"(l - a)^"" . (41) 



apj A/27r«(l - a)p 

It is straightforward to show that the function g satisfies the bound g{a) > 1/2 and 
therefore the previous binomial coefficient grows like g{a)~''^ < 2^. As a consequence, 
the "classical term" dominates the asymptotics in the sense that: 

for a < 1/2. This result shows that the quantum corrections (due to the finiteness 
of k) do not affect the asymptotic expansion of N2{p) neither at the leading neither 
at the subleading order. 



3. Entropy of the SU{2) Black Hole 

Here we adapt the techniques used in [10] and firstly introduced in [8] to compute 
the entropy of a black hole. We are not going into the mathematical details of 
these techniques which has been very well exposed in and which are in fact very 
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well-known in the domain of probabilities and used to understand some properties 
of random walks. We propose to reproduce these results in a more "intuitive" or 
physical way: we will omit many mathematical details which appear in a first time 
non necessary. In particular, we show that it is not necessary to go to complex 
analysis and number theory to get the asymptotic behavior of the entropy. 

Before going to the details of the entropy calculation, let us briefly recall how 
black holes are described in loop quantum gravity and how we compute the entropy. 
In the context of LQG, a local definition of a black hole is introduced through 
the concept of isolated horizons (IH), regarded as a sector of the phase-space of 
GR containing a horizon in equilibrium with the external matter and gravitational 
degrees of freedom. This local definition is used for the black-hole entropy calculation 
since the quantization of such a system allows to define a Hilbert space which is 
the tensor product of a boundary and a bulk terms. The entropy of the IH is 
then computed by the formula S = tr(pj^ logp^ji^), where the density matrix p^jj is 
obtained by tracing over the bulk d.o.f., while restricting to horizon states that are 
compatible with the macroscopic area parameter a. Assuming that there exists at 
least one solution of the bulk constraints for every admissible state on the boundary, 
the entropy is given by S* = log(A^(a)) where N{a) is the number of horizon states. 
Since the theory on the horizon is associated to Chen-Simons theory with punctures, 
the entropy calculation problem boils down to the counting, in the large horizon area 
limit, of the dimension of the Hilbert space ([2]), which, for a given configuration of 
punctures spins d = {di, ■■■ ,dp), is expressed by the formula (l32l) . 

3.1. The Laplace transform method: basic idea 

The Laplace transform method allows in certain cases to obtain the asymptotic 
behavior of a function F[p) for large values of p. For simplicity, we assume that 
the function F is defined for p integers but the method applies in the case where F 
is a function of a real number x. The idea consists first in considering the Laplace 
transform F{s) defined a priori for s > by: 



The Laplace transform appears as a series and therefore might be not defined at 
all or it might be defined for some values of the real positive variable s only. To 
understand if the series is convergent or divergent, one looks at the asymptotic 
behavior of F{p) at large p. To be more concrete, let us propose some examples. 

1. If ~ p" for some real number a, then the series fl43l) is convergent for any 
values of s. 

2. If F{p) ~ 6°"^^ for some positive real number a, then the series fl43p is divergent 
for any values of s and then the Laplace transform is never defined. On 
the contrary, if a is negative, then the series in convergent and the Laplace 
transform is well defined for all values of s. 



oo 




(43) 



p=0 
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3. If F{p) ~ e"^ for some positive real number a, then the series is convergent for 
s > a. In the case where s < a, the series diverges and therefore is ill-defined. 
The case s = a is critical: the convergence properties of F in that regime 
depends on the subleading behavior of F{p). 

We are particularly interested in the last case because in the context of black hole 
entropy the number N{a) of microstates corresponding to a given macroscopic area 
a is exponential in a, namely N{a) ~ e"" and all the problem is to find the coefficient 
a. Let us now come back to the general discussion. If we find a > such that F{s) 
is defined for s > a and undefined for s < a, then we conclude that F{p) ~ e"^ for 
large p. This is more a physical argument than a rigorous proof because we assume 
the asymptotic behavior of F{p). However, it is also possible to prove rigorously 
the asymptotic behavior as it was done in [8] and [10] where the reader can find the 
mathematical details. In the language of probabilities, Sc = a is called the critical 
exponent of F{p). 

Before going further, let us recall that, in the particular case where F{p) ~ e"^, 
it is possible to invert the Laplace transform and to recover to function F{p) from 
F{s) according to: 



where sq > a is a real number. 

In fact, it is possible to extend this technique to obtain also the subleading terms 
in the asymptotics expansion of F{p). To understand this point, we assume that 
F behaves as follows F{p) ~ e°'^Q{p) where Q{p) is an algebraic function whose 
dominant term at large p is Q{j>) ~ p^ where /3 is a real number. We generalize the 
Laplace transform to the following function of the two real variables s and t: 



When t = 0, this function coincides with the standard Laplace transform up to 
the constant -F(O). The point is to evaluate it at the critical value Sc = a. Indeed, 
e~°'PF{p) ~ and therefore the convergence properties of the series F{a, t) depends 
on the asymptotic behavior of p^~*: if /3—t > —1, the series is divergent; if P—t < — 1, 
the series is convergent. Therefore, we proceed as we do to obtain a, we define the 
critical value tc as the minimal value of t such that F{a,t) is well-defined. Thus, 
the coefficient (3 is fixed by /3 = tc — I- 

Of course, we can repeat this technique to obtain recursively all the corrections 
to the leading term in the asymptotic expansion of F{p). But, for this method to 
work, we must know the form of the asymptotic behavior of the function F{p). Once 
we know that the function F{p) ~ e^^p^ at large p, our method allows to obtain the 
critical exponents a and /3. 




(44) 



oo 




(45) 
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3.2. A simple application of the Laplace transform method 

Let us show that we can use this very simple technique to obtain the asymptotic 
behavior of the number of states Nk{d) when all the spins are equal. To that aim, 
we introduce the notation F^ip) = Nk{d) with ji = ■ ■ ■ = jp = j and d = 2j + 1: 

where /iA:(6') is a continuous non singular function obtain directly from ( l32l) : 

2 . 2 ^ sin((2r + 
TT sm(Kfc') 

We omit to mention the dependence in k of ^k{d) and of F(i{p) for clarity reasons. 
Furthermore, we know, from random walks arguments, that the asymptotics of -Frf(p) 
is dominated by its classical part only; in other words, we consider r = in the sequel. 
We will discuss the quantum corrections later. Now, we assume that Fd{p) ~ e""^ for 
large values of p. If the asymptotic behavior assumption is true, then the Laplace 
transform F^i^s) = Yl^=o^~^^-Fi^^'P) is well-defined for s > a and not-defined for 
s < a. To obtain the critical exponent a, we need to simplify the expression of 
F(i{s). To do so, we exchange the sum defining and the integral over 9 in the 
definition of Fd (H6l) . We obtain the following expression: 



Fdis) = - d6 sm^6' > e " — 

^Jo t^A sin^ ; 



(4J 



p=0 

Then, assuming that s is large enough, we perform the sum over p and then: 

Fd(s) = - r de sm^eZ(s,e) with Z(s,e) = (i-e-'^^^^^ 
n Jo \ sm e 



(49) 



The next step consists in analyzing the structure of the singularities of Fd{s). It is 
clear that its singularities come from the poles of the function Z{s, 9) viewed as a 
function of 9. We immediately see that Z{s,9) admits a pole (viewed as a function 
of 9) if and only if > d~^, i.e. s < logd. More precisely, we have the following: 

1. When s = logd, then Z{s, 9) admits an unique pole which is 9^ = 0. 

2. When s < logd, then Z{s,9) admits also at least one pole 9o ^ 0[7r]. At the 
vicinity of 6'o, the (inverse of the) function Z behaves as follows: 

Z~^{s, 6^0 + e) ~ — e((icotan((i6'o) — cotan(6'o)). 

As a consequence, the integral (1481) is divergent. 
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As expected, F{d, s) is defined only for s > logrf. Tlierefore, Sc = logd is tlie critical 
exponent and we have the asymptotic behavior: 

Fd{p) ~ dP . (50) 

Before computing the sub-leading corrections, let us make some important remarks. 

In the first remark, we come back to the exchange of the sum over p and the 
integral over 6 in the computation of the Laplace transform fjlHj) . This step can be 
justified in our case but this is not always the case. More precisely, the exchange 
makes sense if the sum over p is defined, namely if \e~'^ sm{d6)/ sm6\ < 1 for all 6. 
This is exactly the condition we obtained to compute the critical exponent. 

In the second and last remark, we come back to the "quantum corrections" of 
Fd{p)- We know that the number of states is given by: 



FM = 




The first term is the classical contribution and the rest are the quantum corrections. 
We want to compute the asymptotic behavior of these corrections at the large p 
limit, and we take, at the same time, k large as well with a fixed ratio p = p/2k. 
As a consequence the number r = [p] remains fixed in this limit. Therefore, the 
contribution of Fd{p) corresponding to g 7^ reads: 




The calculation of the Laplace transform is more involved in that situation due to 
the presence of a fast oscillating function in the integrand. A consequence is that 
the naive exchange of the infinite sum over p and the integral over 6 is not justified. 

Now, we go further and compute the sub-leading terms. For that purpose, we 
concentrate on the classical contribution only (formally we take t = 0) and we define 
the general Laplace transform of F^^p) P^jl : 

00 

F,{s,t) = Y,e-^'p-'F,{p) . (53) 
Permuting the integral with the sum leads to the expression: 

Indeed, we recognize the polylogarithm function Lit(z) defined for any couple of 
complex numbers {t,z) such that |2;| < 1 by the series: 

00 

Li,(^) = J2p-'zP. (55) 
p=i 
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Following the general idea we described above, we first evaluate -Frf(s, t) at the critical 
value Sc = log d. Then, we look for the critical value tc such that -Frf(sc, t) is well- 
defined for t > tc but not defined if t < tc- We know that Frf(sc,t) might be not 
defined because of the singularity of the integrand Li^ ^^^j^ j at 6* = 0. To compute 
tc, we analyze the behavior of the integrand around ^ = 0: 

1 - 

~ Ut{e—' ) 

where r{t) is the Gamma function. We used some asymptotic properties of the 
polylogarithm function. As a consequence, at the vicinity of 6* = 0, the integrand of 
(13^ behaves as: 

r(i - 1) ^^(-1) ~ ^ r(i - 1) (56) 

because fik{0) — "^6"^ for small. Therefore, the integral over 6 is defined when 
2t > —1 and not defined when 2t < —1; then the critical value of t is tc = —1/2. As 
a consequence, we obtain the value of the second critical exponent 

/3 = -1 + = -3/2 

which is independent of the dimension d. Finally, we can establish that: 

3 

Fd{p) ~ dPp~^/^ then log Fd{p) = plogd - - logp + C(l) . (57) 

In particular, we recover from another method the asymptotic behavior of F2{p) 
given in (1361) . This asymptotics has been obtained in [13] from random walks argu- 
ments. The case d = 2 has also been consider in [T^ . 

3.3. Asymptotic of the entropy 

We apply, the method illustrated above to compute the asymptotic behavior of a 
spherically symmetric and a distorted SU{2) black hole entropy in loop quantum 
gravity. 

a) The spherically symmetric Black Hole 

The calculation of the entropy of the spherically symmetric SU{2) black hole has 
been done precisely in [TU] and has been investigated earlier in [12] when the level in 
infinite. In [HI [7] , we have shown that the level k and the area a can be considered 
as independent variables. For that reason, we are going to reformulate the results 



Li, 



sin(d9) 
dsin 9 



19 



obtained in [TO] when k is finite focussing on pliysical arguments and avoiding tlie 
number tlieory and tlie complex analysis aspects. These aspects are not necessary 
to get the main ideas and the main results if one assumes that the number of states 
grows exponentially with the area. 

The entropy S{a) = log N[a) of a spherically symmetric black hole of macroscopic 
(adimensionalized) area a is defined from the number of states 

N(a) = ±'fna- £kV(*ZM±i))^^y) , (58) 

p=0 di,--- ,dp 

The finiteness of the level k appears in two different places: in the sums which run 
from 2 to + 1 and in the expression of Nk{j). It is important to note again that 
we will consider k and a as independent variables and we will study the entropy for 
large a but finite k. This consideration allows us to define the Laplace transform of 
N{a) 

N{s) = / dae-^'Nia) . (59) 



One difference with the previous Section is that now the variable a is continuous. 
After some calculations and assuming that s is large enough, we end up with the 
following expression: 



= E i: r de ^,{6) lf['^^e-^Vi^^^^^ (60) 
= f; r<(9M«)fi:|;fe-4V<3^y (61) 

p=0 \d=2 ) 

= ^^...^(l-t^^^^e-^y' (62) 

Again, we exchanged the sums over the variables and p with the integral over Q. 
We proceed as in the previous section and conclude that the critical value of s is the 
highest zero of the function 

l-y-!i!^il±ll^e-tv/^ (63) 
^ sin 61 

d=\ 

which is reached for = 0. Therefore, the critical exponent a is the unique solution 
of the equation 



l-^(rf + l)e-tV^(^ = 0- (64) 
d=\ 

With that definition, a depends on the level k. For increasing values of the level /c, 
the solutions a of the previous equation reach fast an asymptotic value as plotted 
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Figure 6: Increasing critical values of the exponent a for different A; G N. Already 
for A; > 4 an asymptotic value is reached. 



in Figure [6l The asymptotic value coincides as expected with the value ctoo found in 
PH] when k ^ oo. Furthermore, we can estimate how the difference Aa = Ooo ~ « 
decreases when k increases. Aa decreases exponentially with k in the sense that it 
exits two real positive constants A and B such that 

Ittfc-aool < Ae~^^ 



Next, we compute the sub leading corrections. The difficulty of this problem is 
that we could not find a way to put the generalized Laplace transform in a "suit- 
able" form as it was the case for the "toy-example" we considered above. We will 
nonetheless circumvent this difficulty as follows. First, we will evaluate N{s) at the 
vicinity of the critical value a, i.e. s = a + e for e small, and then we will compute 
the generalized Laplace transform N{a, —t) at the critical value a as follows: 



d'N{a + e) 



k=0 



assuming that t is an integer. Finally, we will see that it makes sense to extend the 
obtained formula to half-integers (and also real numbers in fact) which will allows 
us to extract the sub-leading corrections to the entropy. 

Let us start, as announced, by the following calculation: 

N{a + e)^ rd9^,{e) (l-j^ ^^^^^ ^^^ ^"^^^(^ ' '7;Vd(dTY))] (65) 
Jo \ smfc/ Z J 

We know that the singularity of N{s+e) when e goes to zero is due to the singularity 
of the integrand 
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in (l65|l when 9 goes to zero. Therefore, we concentrate on the behavior of fs{0) at 
the vicinity of = 0: 



j,{Q) ~ (2r + 1)— ( + ^) V^^(^ + 2) + + + 2))e-tV^(^ 

Now, we start from this expression to study the singularities of the generahzed 
Laplace transform N{q., t) evaluated at the critical value. Indeed, when t is a positive 
integer, we can compute: 

~ d'Nja + e) , 

N{a,-t) = — 1^=0 (66) 

The last quantity is expressed as an integral over the variable 6 whose eventual 
singularity is due to the behavior of the integrand around 9 = given by: 



f ELi^(^+l)(^ + 2)e-tV^) 



t+i 



-2t 



We assume that the behavior of the integrand of N{a, —t) remains the same even 
if t is any real number. As a consequence, N{a,t) is singular when 2t > 1 and the 
critical value of t is tc = 1/2. Then, the critical exponent /3 = —tc — 1 = —3/2 and 
we recover the asymptotic expansion [TO] : 



N{a) ~ e"^a-^/'^ for large a, (67) 

with a given as in Figure [61 Again, the finiteness of k does not modify the sub- 
leading corrections when k is large. 

Let us finish this section by two remarks. The first one concerns the effect of 
the finiteness of the level k in the entropy. As we have just seen, k does modify 
the leading but does not modify the subleading corrections of the entropy. In that 
sense, the logarithmic corrections seems to be universal and independent of the 
Immirzi parameter even in the SU{2) spherically symmetric black hole. The second 
one concerns the techniques we used: our calculations have to be viewed more as 
"physical" arguments than rigorous proofs of the asymptotic behavior of the entropy. 
The nice point is that we can recover the "right" results very easily without entering 
too much into technical aspects. 



b) The distorted Black Hole 

In the distorted case [7J, the black hole is described in terms of two commuting 
Chern-Simons theories associated to the same level k. As in the symmetric case, 
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the area a of the black hole and the level k are considered as independent variables. 
At the fundamental level, the description of the distorted black hole in terms of 
microstates is rather different from the description of the symmetric black hole. 
Indeed, each puncture colored with a SU{2) representation j coming from the bulk 
decomposes into two SU{2) representations and j~ when it crosses the black hole. 
A macroscopic state is therefore characterized by the number p of punctures and a 
family {jt, , € of 3p representations of SU{2) such that 

satisfy the triangular inequality for each i and an additional constraint. 

To describe this additional constraint, we associate canonically the SU{2) gen- 
erators J]_, Jl, and J* to each puncture: the Casimir of these operators J|. and 
are fixed by the representations j± and j in the standard way. The constraint reads 

C\p) = T_-J\- a{J\ + Ji) = 0. (68) 
with ^ ^ 

a ^ (69) 

Now C* and = + Ji + J* = (implicitly imposed above) cannot be simul- 
taneously strongly imposed as they do not form a Lie algebra. One has to impose 
them weakly and there are two possibilities. 

In order to see this let us exploit the fact that there is a strict analogy with the 
way the simplicity constraints are imposed in the EPRL-FK model [2T1 120] . Observe 
first that equation (!68|) has the very same form of the linear simplicity constraints 
of the EPRL-FK models where the role of the Immirzi parameter is here played by 
a. 

The first possibility of weak imposition consists of taking 

j± = (l±a)j72 (70) 

implying 

j=j++j_. (71) 

It can be checked that this choice is consistent with alpha as given in ( l69l) . With 
this then one can check that for an admissible state \ip) one has 

which vanishes in the (semiclassical) limit /i — )■ 0, j — )■ oo with hj kept constant. 
Moreover, one has that 

{<P\CM) = (72) 

for arbitrary pairs of admissible states. In other words, in this first possibility the 
constraint C* are satisfied strongly in the semiclassical limit, and weakly in the sense 
of matrix elements in general. 

The second possibility is not to impose the condition (ITUj) and leave j± completely 
free and only constrained by the triangular inequalities with j. In that case it has 
been shown [22] that ( 1721) is still satisfied. This second possibility is still compatible 
with the classical limit but is weaker than the previous one. 
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Entropy calculation A 

In this case we impose condition ( ITUj) and hence j = j+ + j_. Such a state 
characterizes a black hole of macroscopic area 



1 ^ / 



£=1 



in unit of ip. As a consequence, the number of microstates N{a) associated to a 
distorted black hole of area a is given by the formula: 



where df = 2jf + 1. Following the steps of the previous Section, we introduce the 
Laplace transform N{s) of the number of states N{a). It is given by: 

oo k+l 



n J 1 Jo Jo 



p=0 d+,d 
^ P 



p . fd+e+^ . ,d- 



1=1 sin ^ sm ^ 



where the sums run over the families = [di, ■ ■ ■ ,dp) of representations dimen- 
sions. Following the same strategy as in the spherically symmetric case, the previous 
expression can be simplified as follows: 



N{s) = J2 de^i^kie^) / de~ 

p=o Jo Jo 




— e 2 



' - ' '>^(d++d--i){d++d-+i) 



sm — 



with 



In the definition of D^, the sums run over e [0, A;], which are related to the spin 
variables by the relation = due to the changing of variables. 

As in the spherically symmetric case, we conclude that the critical value of s is 
the highest value for which D^, viewed as a function of the angles 9"^, admits a zero. 



24 



It is reached when admits one zero at 9"^ = 0. Therefore, the critical exponent 
a is the unique solution of the equation 

k 

1 _ ^ + + i)e-fV(4+^7+i)(4+rf7+3) = . (75) 

d±=0 

The exponent a depends on k as in the spherically symmetric case and the numerical 
solution has been plotted in Figure [3 



1.5 ^ 
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Figure 7: In the figure, we have plotted the values of the exponent a as function of 
A; G N for the first integers. The plot shows as, similarly to the spherically symmetric 
black hole counting, in the case of a weak imposition of the constraint C" = through 
the relation j = j+ + j_, an asymptotic value for a is quickly reached as A; > 4. 

Entropy calculation B 

Let us now concentrate on the case of the weaker imposition of the constraint 
C* = 0, where all pairs of admissible states are taken into account. In this case, the 
black hole of macroscopic area is 

I P 

a = - Y^^J{dt-l){dt + l) 

£=1 

in unit of £p. As a consequence, the number of microstates N{a) associated to a 
distorted black hole of area a is given by the formula: 

p=0 di,--- ,dp 

k+1 / p \ _ _ 

E [U^iji,jl,Ji)] Nkij^)N,{r), (76) 

d±...,d±=l \^=1 / 

where df = 2jf + 1 and we recall that, in order to implement the admissibility 
condition, Y{je, j'l , j^) = 1 if satisfy the triangular inequality, it vanishes 

otherwise. 
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Following the steps of the previous case A, we introduce the Laplace transform 
N{s) of the number of states N{a). It is given by: 



n J J- "'0 JO 



p=0 d d+,d 



sm \ sm 



sm ^ sm ^ 



Following the same strategy as in the case A, the previous expression can be sim- 
plified as follows: 



sm( ^-^) sm' " 



E E ^(^■^'^■Z'- 



2 > 



with now 



J j+ J 1 sm ^ sm -TT 

d d+.d-=l 2 2 



sinf (^^^) sinf (^^^ 



Dfr,r;.) = l-EE^(^-^-^-^^7) \ ' . I ^ -tv^.(77) 

Again, the spins variables j"^ are related to the sums variables by = due to 
the changing of variables. The variable d is also related to the spin variable j by 
d = 2j. 

Similarly to the previous cases, we conclude that the critical value of s is the 
highest value for which Dj^, viewed as a function of the angles 6"^, admits a zero. It 
is reached when admits one zero at 6^ = 0. Therefore, the critical exponent a 
is the unique solution of the equation 



1 - E E id-" + l){d- + l)e-fV^ = 0. (78) 

d d±=0 

The exponent a again depends on k and the numerical solutions has been plotted 
in Figure El 

To have a more physical intuition of the behavior of a as a function of k, let us 
assume that all the spins j are fixed to 1/2: this means that the edges of the spin- 
network in the bulk which intersect the black hole surface are colored by 1/2-spins. 
This assumption will give us the "shape" of the function a{k) for large values of k 
since the main contributions to the entropy come from small values of the bulk spin 
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Figure 8: In the plot, the circles indicate the values of the exponent a as function 
of A; G N for the first integers; the squares represent values of the function clog/c for 
the same values of /c, where c = 2\/3 + 

jg. Let us call ai/2 the value of a where only j = 1/2 spins contribute and a;i/2 
satisfies: 

^ ^ (rf+ + l)(ti- + l)e-"V2^ = 1. (79) 
A straightforward calculation leads to the following expression relating ai/2 and k: 

n=l 

As a consequence, in the limit where k is large 

ai/2 ~ 2v^ log k 

which means that a grows logarithmically with k. Further evidence of this behavior 
of a is given by the numerical solution of eq. (!75|) plotted above (Figure [H]) where all 
spins are taken into account and not only 1/2 spins. Note that the numerical value 
c in the asymptotic formula a. ~ clog/c is such that c = 2-\/3 + OiXjk) as expected. 

Let us now concentrate on the sub-leading corrections. We proceed exactly as in 
the spherically symmetric case: we first evaluate iV(s) at the vicinity of the critical 
point a, i.e. s = a + £ for a small e\ we can then study the singularities of the 
generalized Laplace transform A^(a, t) evaluated at the critical value through the 
relation ( l66i) . by expanding the integrand fs{9~^,9^) around 6*+ = = 9^; finally, 
we look at the maximal value of for which N{a,t) is well-defined. The critical 
exponent /3 is then given by /3 = —tc — 1. More precisely, in the distorted case we 
have 



de' (^+2(*+i) + ^-2(t+i)) p2{t+i) (sin2(*+i)(<^) + cos2{*+i)((^)) ' 
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where in the last equahty we have changed to polar coordinates 9^ = psin{(f), 
9^ = pcos{ip). From the previous equation we get that 



N{a,-t) 




Consequently, N{a,t) is singular when t > 2, in this case the critical value of t is 
tc = 2. Therefore, the critical exponent now is /3 = —tc — 1 = —3 and the asymptotic 
expansion reads: 

N{a) ~ e'^'a-^ , (81) 

where a is given in Figure [7]or|8] according to the prescription that defines the allowed 
states. From the previous expression for the Laplace transform, we see that, in the 
distorted case, the constant factor in front of the logarithmic corrections becomes 3. 

Conclusion 

This paper has been devoted to the calculation of leading and sub-leading terms 
of the SU{2) black hole entropy in Loop Quantum Gravity when the black hole 
is spherically symmetric [HI E] and when it is distorted [7J . To reach this aim, we 
derived first, by means of the recoupling theory of the quantum group Ug{su{2)), a 
new integral formula, resulting to be very useful, for the dimension of the Hilbert 
space of SU{2) Chern-Simons theory on a punctured two-sphere, which enters the 
definition of the Hilbert space of the SU{2) spherically symmetric and distorted 
black hole as derived in [HI [7] . Successively, we revised the technique of the Laplace 
transform method, exposed in detail in [10] and firstly introduced in |8j, to study 
the asymptotic behavior (in the large area limit) of the entropy associated to these 
two statistic mechanical ensembles. 

The entropy of a SU{2) spherically symmetric black hole has been already studied 
in [in] when the level k is infinite. Here, following a paradigm-shift introduced in 
[7], we considered the case where the level k of the Chern-Simons theory and the 
macroscopic area of the black hole a are independent variables and we studied the 
effect of a finite k. We showed that, if one takes into account the finiteness of 
the level, the entropy of type I isolated horizons is not modified at least up to the 
subleading corrections, therefore recovering the results of [TO]. Moreover, the critical 
exponent of the leading order a, which is now a function of the level k, reaches fast 
an asymptotic value for large k, as shown in the plot in Figure 6. 

Concerning the entropy of a distorted SU (2) black hole, this is something which 
has never been studied before. In this case, for each puncture coming from the bulk, 
there are two punctures associated to it on the horizon and the Hilbert space becomes 
now the direct product of two SU (2) Chern-Simons Hilbert spaces with same level k 
[7]. The SU{2) symmetry is implemented by the insertion of an intertwiner between 
the three punctures (one from the bulk and two from the horizon), therefore, the 
area constraint still plays an important role. Using the techniques developed for the 
spherically symmetric case, we have performed the counting of the enlarged Hilbert 
space number of states and shown that the entropy is again proportional to the 
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horizon area to the leading order. In the distorted case, one can distinguish two 
different models according to the way second class constraints are imposed weakly. 
In the strongest imposition of the constraints the results do not differ in a qualitative 
sense from those obtained in the spherically symmetric case. However, if the second 
class constraints are only imposed weakly, in the Gupta-Bleurer sense, then the 
critical exponent a does not go to an asymptotic value for increasing values of the 
level k but grows logarithmically with it, as shown in Figure [HI In that sense, our 
model is consistent with "any" value of the Immirzi parameter. 



Acknowledgments 

This project was partially supported by the ANR. D.P. was supported by the Marie 
Curie EU-NCG network. 



References 

[1] A. Ashtekar and J. Lewandowski, Class. Quant. Grav. 21, R53 (2004) 
[arXiv:gr-qc/040 4018| . 

C. Rovelli, "Quantum Gravity," Cambridge, UK: Univ. Pr. (2004) 4 55 p. 

T. Thiemann, "Modern canonical quantum general relativity," Cambridge, 
UK: Umv. Pr. (2007) 819 p. 

A. Perez, "Introduction to loop quantum gravity and spin foams," 
[a?Xiv:gr-qc/040906l] . 

[2] 
[3] 



[4; 

[5 

[6 

[7: 



A. Ashtekar and B. Krishnan, "Isolated and dynamical horizons and their 
applications," Living Rev. Rel. 7 (2004) 10 | arXiv:gr-qc/0407042| . 

A. Ashtekar, J. C. Baez and K. Krasnov, Adv. Theor. Math. Phys. 4 (2000) 1 
|arXiv:gr-qc/0005126] . 

C. Beetle and J. Engle, "Generic isolated horizons in loop quantum gravity," 
Class. Quant. Grav. 27 (2010) 235024 [ar Xiv: 1007.2768 [gr-qc]]. 

J. Engle, A. Perez and K. Noui, "Black hole entropy and SU(2) Chern-Simons 
theory," Phys. Rev. Lett. 105 (2010) 031302 [arXiv:0905J168j [gr-qc]]. 

J. Engle, K. Noui, A. Perez and D. Pranzetti, "Black hole entropy from an 
SU(2)-invariant formulation of Type I isolated horizons," Phys. Rev. D 82 
(2010) 044050 [ arXiv:1006.063 4 [gr-qc]]. 

A. Perez and D. Pranzetti, "Static isolated horizons: SU(2) invariant phase 
space, quantization, and black hole entropy," ar XivilOll. 29611 [gr-qc]. 

K. A. Meissner, "Black hole entropy in loop quantum gravity," Class. Quant. 
Grav. 21 (2004) 5245 [arXiv:gr-qc/0407052|. 



29 



M. Domagala and J. Lewandowski, "Black hole entropy from quantum geom- 



etry," Class. Quant. Grav. 21 (2004) 5233 | arXiv :gr-qc/040705l|. 

[9] A. Ghosh and P. Mitra, "Fine-grained state counting for black holes in loop 
quantum gravity," Phys. Rev. Lett. 102 (2009) 141302 [arXiv:0809.4T70l [gr- 
qc]]. 

A. Ghosh and P. Mitra, "Counting black hole microscopic states in loop 
quantum gravity," Phys. Rev. D 74 (2006) 064026 [arXiv:hep-th /0605125|. 

A. Ghosh and P. Mitra, "Counting of black hole microstates," Indian J. Phys. 
80 (2006) 867 |arXiv:g r-qc/0603029] . 

A. Ghosh and P. Mitra, "An improved lower bound on black hole en- 
tropy in the quantum geometry approach," Phys. Lett. B 616 (2005) 114 
[arXiv:gr-qc/04110 35] . 

A. Ghosh and P. Mitra, "A bound on the log correction to the black hole area 



law," Phys. Rev. D 71 (2005) 027502 arXiv:gr-qc/0401070|. 



[10] 1. Agullo, J. Fernando Barbero, E. F. Borja, J. Diaz-Polo and E. J. S. Vil- 
lasenor, "Detailed black hole state counting in loop quantum gravity," Phys. 
Rev. D 82 (2010) 084029 |arXiv:1101 .3660 [gr-qc]]. 

I. Agullo, G. J. Fernando Barbero, E. F. Borja, J. Diaz-Polo and E. J. S. Vil- 
lasenor, "The combinatorics of the SU(2) black hole entropy in loop quantum 
gravity," Phys. Rev. D 80 (2009) 084006 ^arXiv:0906.4529i [gr-qc]]. 

J. F. Barbero G. and E. J. S. Villasenor, "On the computation of black 
hole entropy in loop quantum gravity," Class. Quant. Grav. 26 (2009) 035017 
[arXiv:08 10.1599 [gr-qc] ] . 

J. F. Barbero G. and E. J. S. Villasenor, "Generating functions for black 
hole entropy in Loop Quantum Gravity," Phys. Rev. D 77 (2008) 121502 
[arXiv:0804.4784l [gr-qc] ] . 

I. Agullo, J. F. Barbero G., J. Diaz-Polo, E. Fernandez-Borja and E. J. S. Vil- 
lasenor, "Black hole state counting in LQG: A number theoretical approach," 
Phys. Rev. Lett. 100 (2008) 211301 [ arXiv:0802.4077l [gr-qc]]. 

[11] H. Sahlmann, "Entropy calculation for a toy black hole," Class. Quant. Grav. 
25 (2008) 055004 jarXiv:0709.0076l [gr-qc]]. 

H. Sahlmann, "Toward explaining black hole entropy quantization in loop 
quantum gravity," Phys. Rev. D 76 (2007) 104050 [a rXiv:0709.2433 [gr-qc]]. 

[12] R. K. Kaul and P. Majumdar, "Quantum black hole entropy," Phys. Lett. B 
439 (1998) 267 |arXiv:gr-qc/9801080|. 

R. K. Kaul and P. Majumdar, "Logarithmic correction to the Bekenstein- 
Hawking entropy," Phys. Rev. Lett. 84 (2000) 5255 [ar Xiv:gr-qc/000204 0|. 

S. Das, R. K. Kaul and P. Majumdar, "A new holographic entropy bound from 
quantum geometry," Phys. Rev. D 63 (2001) 044019 |ar Xiv:hep-th/000621l] . 



30 



[13] E. R. Livine and D. R. Terno, "Quantum black holes: Entropy and entangle- 
ment on the horizon," Nucl. Phys. B 741 (2006) 131 |arXi v:gr-qc7b 508085|. 

L. Freidel and E. R. Livine, "The Fine Structure of SU(2) Intertwiners from 
U(N) Representations," J. Math. Phys. 51 (2010) 082502 |arXiv:0911.3'553l 
[gr-qc]]. 

E. R. Livine and D. R. Terno, "The entropic boundary law in BP theory," 
Nucl. Phys. B 806 (2009) 715 [arXiv:0805 .2536 [gr-qc]]. 

E. R. Livine and D. R. Terno, "Bulk Entropy in Loop Quantum Gravity," 
Nucl. Phys. B 794 (2008) 138 [arXiv:0706.0985 [gr-qc]]. 

[14] Achucarro A., Townsend P., "A Chern-Simons action for three-dimensional 
anti-de Sitter supergravity theories", Phys. Lett. B 180, 85-100 (1986). 

[15] Witten E., "2+1 dimensional gravity as an exactly soluble system", Nucl. 
Phys. B 311, 46-78 (1988). 

[16] Witten E., "Quantum field theory and the Jones polynomial" , Commun. Math. 
Phys. 121, 351 (1989). 

[17] C. Meusburger and K. Noui, "Combinatorial quantisation of the Euclidean 
torus universe," Nucl. Phys. B 841, 463 (2010) |arXiv: 1007.4615) [gr-qc]]. 

[18] Chary V., Pressley A., "A guide to quantum groups", Cambridge University 
Press (1994). 

[19] Verlinde E., "Fusion Rules and Modular Transformations in 2D Conformal 
Field Theory", Nucl.Phys. B 300 3Q0 (1988). 

[20] J. Engle, E. Livine, R. Pereira and C. Rovelli, "LQG vertex with finite Immirzi 
parameter," Nucl. Phys. B 799 (2008) 136 [a rXiv:0711.0146 [gr-qc]]. 

[21] L. Freidel and K. Krasnov, "A New Spin Foam Model for 4d Gravity," Class. 
Quant. Grav. 25 (2008) 125018 [arXiv:0708.1595l [gr-qc]]. 

[22] Y. Ding, M. Han and C. Rovelli, "Generalized Spinfoams"; 
(aBCiv:gr-qc/1011214p ]. 



31 



